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1 Introduction 

In high energy particle colUsions, the interaction volume is often very small 
and large deviations can occur from the thermodynamic limit depending on 
the beam energy. In some cases it may be easy to produce a pion since one 
needs only 140 MeV for its rest mass, however, to produce an anti-proton one 
needs at least twice the rest mass, 1.88 GeV and not simply 0.94 GeV, since 
they can only be produced in pairs. The probability to produce anti-protons 
therefore cannot simply follow the same law (e.g. a Boltzmann distribution) 
as the production of pions. Fortunately, statistical mechanics provides us with 
the tools to take into account constraints like baryon number or strangeness 
conservation. This will be presented in these lecture notes. 

Basically, if the system is large and hot, the corrections are negligible. 
In Pb — Ph collisions at the CERN/SPS collider they are negligible because 
the system is large enough but for p — p collisions at the same energy these 
corrections are large and must be taken into account because the final system 
is too small. 

The main statistical concepts are collected in Ch. 3. In chapter 3, a sta- 
tistical method for taking the exact strangeness conservation into account is 
presented. This is applied in Ch. 3.1 to describe the particle production as 
observed at the GSI in Ni-|-Ni collisions. The exact treatment of quantum 
numbers is extended to include the strangeness, baryon number and charge 
in Ch. 4. A comparison of numerical results with the AGS E802 data is given 
in Ch. 4.1. Finally in Ch. 5, the generalization of the method to arbitrary 
internal symmetry is reviewed. 



2 Quantum Statistical Concepts in Brief 

Throughout this paper we follow the usual convention of natural units, so the 
speed of light, Planck constant and Boltzmann constant have values c = 1, 
h=l and fc = 1, respectively. 
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All the physical information about a collection of particles is contained in a 
density operator, p. The average of an observable A in this statistical ensemble 
is calculated as {A) = ix{pA), where A is a Hermitian operator corresponding 
to the observable. Using the density operator, one defines the entropy of the 
system considered as 5' = — tr(/51n/3). In any system in nature, the entropy 
is known to tend to its maximum, so one has to find a representation of the 
density operator satisfying this condition. In thermodynamical equilibrium, 
the average occupations of different quantum states do not change in time, 
so dp/dt = 0. The density operator satisfies the equation of motion of the 
form idp/dt = — [p, H], where H is the Hamiltonian of the system. Thus, the 
thermodynamical, stationary density operator is diagonal in the basis formed 
by Hamiltonian eigenstates. 

The choice of constraints used in maximizing the entropy defines the type 
of the statistical ensemble obtained. The closed system with fixed energy, 
E, volume, V, and number of particles i, Ni, is a microcanonical ensemble. 
System in heat bath (the ensemble average of energy {E) = tr(/5-ff), V and 
Ni are conserved) corresponds to canonical ensemble. Further, if we let the 
particle number Ni fiuctuate such that the average tr(/57Vi) is conserved, we 
obtain a grand canonical ensemble. 

Maximization of entropy using the canonical boundary conditions leads 
to the density operator of the form p = e^^^ /Z, where /3 is the inverse of 
temperature T, and Z is the canonical partition function, 

Z = tre-'5* = ^e-''^^W. (1) 

i 

Here i labels the different quantum states in the system and Ei{N) is the 
eigenvalue of the A'' particle Hamiltonian. In the last step the trace is ex- 
pressed in the basis of Hamiltonian eigenstates. Once knowing the correct 
partition function, one is able to calculate the thermodynamical quantities 
describing the system. For example, the average energy is {E) = T'^d In Z/ dT. 
Choosing the grand canonical boundary conditions yields 

PG = e- IZg- 
Here we have employed the grand canonical partition function, 

Zg{T, {AO, V) = tre-^^"-^'^'^ = [] E (2) 

i Ni 

where Aj = e'^'** is the fugacity of the particle species i and Z^^ is the 
particle canonical partition function of the species i. Chemical potentials /Zj 
take care of particle number conservation in an average sense. In this work, 
we are mainly interested in the mean particle numbers in the grand canonical 
system, so we employ frequently the equation 

m^K'-^. (3) 
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In properly quantized, finite volume system the partition function is often 
rather difficult to compute. In this work we never know the exact geometry 
of the system, so we settle for a finite volume, V, sample of the infinite vol- 
ume system. Thus, the summation over discrete quantum states in partition 
function is changed to simple phase space integration over continuum. The 
one particle canonical partition function of particle i is now given by 



(4) 



where gi is the spin degeneration factor, and rui is the mass of the particle i. 
Using the previous result and taking care of the correct ocupation of quantum 
states, the grand canonical partition function can be written in the form 



\nZGiT,{X,},V)=J2! 



V 
'2^ 



dpp^ In 1 + Ty.A.e-^'VP^ 



(5) 



where rj is the statistics factor: 77,; = 1 for fermions and 77^ = —1 for bosons. 
Now we can write the mean particle number as 



V 
2^ 



dpp 



(6) 



In a rare gas (i.e. Boltzmann) limit, which is mostly applied here, we just put 
the statistics factors 77^ = in particle numbers formula, or let the possible 
occupation of one particle states be only one to obtain 



InZc = ^A,Zi. 



(7) 



In the relativistic, multispecies gas, where the conservation of number of 
distinct particles is not the main interest, we associate the chemical potentials 
to conserved quantum numbers. Given that baryon number B, strangeness 
S and electric charge Q are conserved averagely in a relativistic hadron gas, 
the particle numbers (p) are 



V 

2^ 



dpp^ 



+ Vi 



(8) 



where Bi, Si and Qi are the quantum numbers of individual particle species. 
The average net quantum number, say baryon number, is a sum over particle 
numbers weighted by the chosen quantum number. 



i 



(9) 



4 J. Cleymans et al. 



3 Exact Strangeness Conservation 



Let us consider first a gas composed of neutral kaons and antikaons and 
request that the overaU strangeness be zero. The canonical (with respect to 
strangeness) partition function is given by 



Zs=o — 



OO / OC 

n 



(10) 



where et is the energy and rii the number of if°'s in the state denoted by i. 
The Kronecker delta ensures that the overall strangeness is zero: n^o — n-^^o = 
X]fco('^* ^ ^i) ~ 0. By including the l/riil and 1/ni! one gets the sums 
over all distinct states. Replacing the sums over single particle levels by the 
Boltzmannian momentum integrals and using the Dirac representation 



S{n — m) = 



27r 



■In 



dae*("- 



(11) 







of delta function, one obtains the following result 

1 f^^ 
Zs=o ^ -z- da exp 



X exp 



V 
V 

2^ 



dpp e 



dpp e 



(12) 



where ej^o — -^/p^ + . Applying the notation of single particle partition 
function (0) gives 



1 r^'' / 

Zs^o ^^J^ da exp [z^o e"" + ZL, t 



(13) 



By expanding the exponentials in power series it is easy to perform the a 
-integration to obtain 



This series converges to a modified Bessel function 



Zs=o — -^o(2y Z^oZl-o) 
which is generally defined by 

OO 



1=0 



II {l + iyy. \2 



(14) 



(15) 



(16) 
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The generalization of the calculation to a gas containing any particles i 
carrying strangeness Si = 0, ±1 is straightforward. For this case one obtains 

Zs=o^^ f dacxp(Me'"+AA_ie-*"+AAo), (17) 
'^''^ Jo 

where Af^ stands for the sum of single particle partition functions of particles 
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+ zL, + ... 
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(18) 
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f . . . + 
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In this case, where the strangeness is being treated canonically and the baryon 
number and electric charge grand canonically, the Boltzmannian one species 
partition function is 



Z. 



B ^'Q 



27r^ Jo 



(19) 



Performing the power expansion of the exponentials and the integration in 
eq. (H^) we are left with 



Zo -^i^iM-if = ZoIoi2y^N\Mi;), 

p=0 



1 



(20) 



where Zq is the grand canonical partition function for particles having strangeness 
zero. To calculate the average abundancies of particles, we substitute the fic- 
titious strangeness fugacity and derive 



which gives 

= Z, 



d In Zs=c 



^X^ 



I M \ /s,(2VVVVVU) 



(21) 



(22) 



Each term in the sum of eq. ( pO[ ) is the product of a strangeness plus one 
and a strangeness minus one particle and one sees the exact strangeness 
conservation explicitly at work. Due to the strict conservation, the number 
of strange particles increases nonlinearly with volume, which is illustrated in 
Fig. (0). 

The method used above and the expression of the partition function Zs=o 
indicate that to impose the strict strangeness conservation, one projects the 
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1.2 



h(X) 




X = 2y'MiM-i (X r 

Fig. 1. Nonlinear volume dependence of strange particle abundances in canonical 
treatment (see eq. pit) compared to grand canonical case {X = 1). 



grand canonical partition fmiction Zg{T, Xb, Xq, Xs) onto the state with 
strangeness S, 



= ^ dae-'^"ZG(r, Xb,Xq, Xs), (23) 

where the fugacity factor As has been replaced by 

As = e^". (24) 

The partition function for a gas containing particles with strangeness 
Si = 0, ±1, ±2 is given by 

Zs^— I dae-^-^^expTMe^" +AA_ie-'" (25) 
27r J a 

+AA3o(ABe"''" + As^e'^a) + A^^- (AsAg^e-''^" + As^Age^^")] , 

where the sums A/h-o and Ms- include also the heavier resonances carrying 
the same quantum numbers, as the A/±i do. Using the generating function 
for modified Bessel functions Ii, defined by 



exp||U+-j U ^'^iP)^'' (26) 
and expanding again the exponentials in power series we obtain 

Zs^^ r dae-'^'^ V /„(2AAso)AS'e-^2ma (3?) 
zvr /n — 
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OO 

X 



OO OO ^ OO 



p! ^ — ' g! 

n— — OO p—0 q—0 

Carrying out the integration and rearranging the summations the result can 
be expressed in terms of -functions: 

OO OO 

Zs = Zo I^n{2Uso)X^ ^(2A/h-)A'^Aq" 



m— — OO n— — OO 



OO ^ 



pi (-S +p-2m-2n)\ 

p—0 ^ ^ / 



~S-\'2p-2rn-2n / J\f^ \ 2 



oc 00 



m— — oo n— — 00 

5+2m+2n 



+2„(2/^^TArr). (28) 



Using the same techiques the result can be generalized to the case, where the 
il -like hadrons {Si — ±3) are included as well: 

oc oc 

Zs = Zo Y Im{2UEo)yS Y ^n(2A/s-)ASAQ" 



rn— — OO n— — OO 



Y Ii{2Nn-)\'B^Q' 



X 

/ — — OO 



I M 



Is+2m+2n+3l {2^NiN-i). (29) 



The mean number of hadrons i with strangeness Si in the system becomes 

The modified Bessel functions decrease quickly with increasing absolute value 
of their indecies, so the numerical evaluation of mean particle numbers is not 
cumbersome. 



3.1 Application of Zs 

In a recent paper |^ we have analysed the particle production in GSI SIS 
Ni-|-Ni experiments. We addressed especially the abundance of kaons who 
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can not be described by hadronic gas model in its standard form. Although 
the size of the Ni system is relatively large the corrections due to exact 
strangeness conservation turned out to be crucial at low temperatures, T < 
100 MeV, involved. At these temperatures the width of resonances had to be 
taken into account. A summary of our results for particle ratios is presented 
in table 1. 



Table 1. Particle ratios given by present model compared to experimental results. 
The best fit value, /is — 0.72 GeV, for the baryon chemical potential is used. 



Ratio 



Model 

i? = 4.2 fm i? = 3 fm 

T = 65 MeV 75 MeV T = 65 MeV 75 MeV 



Data 



ratio 



ref. 



K+/K- 


25.7 


22.4 


23.9 


21.1 


K+/-R+ 


0.0071 


0.0339 


0.0027 


0.0132 


cP/K- 


0.103 


0.082 


0.276 


0.212 


n'^ /n 


0.893 


0.895 


0.894 


0.898 




0.008 


0.015 


0.008 


0.015 




0.225 


0.247 


0.224 


0.246 


tt'VB 


0.104 


0.108 


0.104 


0.107 


d/p 


0.129 


0.188 


0.129 


0.188 



21 ±9 

0.0074 ± 0.0021 

0.1 

0.89 

0.037 ± 0.002 
0.195 ±0.020 
0.125 ±0.007 
0.26 



2, 

11 

8| 

5| 

m 

El 



The measured hadronic ratios with corresponding errorbars are described 
as bands in the (T, (Xb) plane as shown in Fig. 1. The intervals of temperature 
and of chemical potential 

T = 70 ± lOMeV 
fiB = 720±30MeV 

give a good fit to the data. The freeze-out radius of i? ~ 4 fm was extracted 
from the volume dependence of the ratios K~^/tt^ and (j)/K~ . 



4 Exact Baryon, Charge and Strangeness Conservation 

In the case of three exactly conserved, additive quantum numbers we start 
from the single particle partition function of particle i, 

V 

^^^J29^^l '^PP^^'^"'^b,,bM,,qM,s.- (31) 

Making use of the integral representation for 5 -functions and the overall 
conservation constraints 

Y,B,=B, -g, ;^5,; = 5, (32) 
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0.02 I 1 1 1 1 1 

0.6 0.7 0.8 0.9 1 1.1 

I^B [GeV] 

Fig. 2. Curves in the {ij,b,T) plane corresponding to the hadronic ratios indicated. 
The interaction volume corresponds to the radius of 4.2 fm, and the isospin asym- 
metry is B/2Q = 1.04. 



the resulting intergral corresponds to a projection of the grand canonical 
partition function onto the desired values of B, Q and S: 

1 p27T p27r p-Itt 

xZg{T,\b,Xq,Xs)- (33) 

Here the fugacity factors have been replaced by 

XB = e'^, XQ = e"^, Xs = e'". (34) 

As the contributions always come pairwisc for particles and antiparticlcs, the 
fugacity factors will give rise to the cosine of the corresponding angle. It is 
useful to group the particles according to their quantum numbers. Leaving 
out charm, bottom and Si > 2 particles wc arc left with ten categories which 
will be labeled by their particle content. For instance, A/fo stands for the 
sum of one particle partition functions of neutral strange and antistrange 
mesons (if", K , K*^, K* , ...), Mk" stands for charged strange mesons 
{K^, K~ , K*~^, K*~ , ...), and Ma stands for all neutral hyperons and 
antihyperons. With these notations the partition function can be rewritten 
in the form 

Zb,q,s = r^^e-'""^ rdcPe-'Q^ /'"dae"^^" (35) 

[■^T^r Jo Jo Jo 
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X exp{2Afn cos ip + 2A/'^c cos (j) + 2AfK0 cos a + iMK" cos((/) + a) 
+2Np cos{tp + 0) + 2J\f^- cos{i}j -(p) + 2Na++ cos{tp + 2(j)) 
+2Na cos{ip - a) 

+2Ne+ cos{\p + (/) - a) + 2Ns- cos{\p - (j) - a)} . 

The integration above can not be done directly due to cosine terms of mul- 
tiple angles. To circumvent this difficulty, we introduce a new angle whenever 
more than one appears. For example, in the term involving Mp we introduce 
an intermediate angle ^ in the following way 

1 = r Ai5{^ + <t>-0= E ^ f dee*^('^+*-«) . (36) 
The application of the integral representation of the modified Bessel function, 

/„(z) = i/ do; 6^'=°''" cos nw, (37) 
allows one to write the partition function in the form 

Zb,qAt.v) = zJw ) 

n^, — — oo/ 

-^-Q+ni+ris-ns+ns-ne+ZnTl^MrO (38) 
X -^-S+ni-n4-n5-n6(2A/i<-o) 
X (2A/'kc ) J„, (2A/'p)/„3 {2Na- ) 
X {2NA)In, {2AfE+)Ine (2AA^- )/„, {2J^A++)- 



The differentiation of cq. (|3^) for particle abundances yields the result 

Zb,q,s 

The evaluation of the canonical partition function with three simultane- 
ously conserved quantum numbers becomes numerically very time consuming 
for large values of B. So far, for systems with i? > 20 we have been forced to 
resort to the grand canonical treatment. 



4.1 Application of Zb,q,s 

In order to compare our numerical results with the E802 experimental data 
shown in Table ^, we estimate the baryon number and charge of the exper- 
imental system via geometrical considerations. Letting Rp and Rt be the 
radius of a projectile and target nucleus respectively, we assume that the 
radii are directly proportional to the cubic roots of the mass numbers Ap 
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and At of the interacting nuclei. In the case of central coUisions the interac- 
tion region is taken to be a cyhnder of radius i?p, length 2y^R^ — i?p plus, 
two remaining spherical segments at the ends of the cylinder with Rt as 
the radius. We find the total number of participating nucleons (or, baryon 
number, B) and the total charge Q to be 



Ap 
At 



(40) 



Q = Zp + Zt {l- 



1 - 



Ap 
At 



(41) 



respectively. 

The comparison of our numerical results with the E802 data (Table 2) 
of the K'^/tt^ and /n^ ratios as functions of the size of the system or, 
equivalently, B are shown in Fig. 3. In the upper figure we observe that the 
theoretical curve for B/2Q = 5/6 approximates the K^/tt^ data best. The 
theoretical curves lie systematically above the data but drop closer as B/2Q 
decreases towards the collision value. The effect of the isospin asymmetry of 
the system is seen also in the /tt^ data comparison. As the ratio B/2Q 
approaches the collision value the theoretical curves begin to approximate 
the data more closely. 



Table 2. Experimental results reported by the E802 collaboration. B and Q are 
calculated using equations M and p. 



Collision 


K+/n+ 


Ref. 


K-/TT- 


Ref. 


B 


Q 


p + "Bee, 

P + ^^Al27 
P + ^""CUM 
P + ™AU197 
"5i28 + ''^Ui97 


7.8±0.4% 
9.9±0.5% 
10.8±0.6% 
12.5±0.6% 
18.2±0.9% 
19.2±3% 


(l 
(l 


15 
~ 

[J 

11 


1 

1 
1 


2.0±0.2% 

2.8±0.3% 
3.2±0.3% 
3.6±0.8% 


1 

1 
1 
1 


w 

131 


3.9 
5.4 
6.9 
9.7 
102.7 


2.3 
3.1 
3.7 
4.5 
44.0 



5 Generalization of the Projection Method 

In this section, we review the projection method generalized to arbitrary 
internal symmetry of the system in addition to C/(l) of strangeness and U{l)x 
U{1) X U{1) of baryon number, electric charge and strangeness. For complete 
derivation, see the original texts of Turko, Redlich and Hagedorn iQ. The 
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Fig. 3. Thermal model expectations for the production ratios /■k'^ and K~ /-k' 
at a temperature of 100 MeV and a baryon density of 0.04 fm~^ compared to 
experimental results from the Brookhaven AGS. The experimental ratios from Si — 
Au collisions {B ~ 103) is moved to B = 21 for the sake of convenience. 



general method is suitable for non-abelian symmetries, such as SU{2) of 
isospin jl^ or angular momentum [|l6j, and SU{5) of color ijl^ as well. 

If the system is exactly symmetric under the operations of internal sym- 
metry group G, the corresponding group generators Qk have the same eigen- 
states as the Hamiltonian. Thus 



[H,Qk]=0, k = l,...,n, 



(42) 
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where n is the number of parameters in the group. Let us define the generating 
function Z as Z = tr[J7(g)e^^^], where U{g) is an unitary representation of 
the group. With the aid of irreducible presentations of [/, this decomposes to 



r(g) 



z^. 



(43) 



Here we used a character x'^{g) and the dimension d{v) of an irreducible 
presentation U'^{g), and the corresponding canonical partition function Z^. 
Using the orthogonality of characters, 



Ma)r{9)r is) = ^u.' 



(44) 



we may compute the canonical partition function once we know the generat- 
ing function: 



Z,^d{u) J dfi{g)x''i9)Z. 
Further investigation of the generating function reveals that 



(45) 



Z = trexp -pH + i ^ Qkjk 
\ fc=i / 



nnE-P 

j=l p=l n 



fc=i 



(46) 



In the last step, we have expressed the trace in the basis of n -particle Hamil- 
tonian eigenstates. The g^''^ are the conserved charges, and the 7^ are the 
variables of the Cartan subgroup of the group G of rank r. Eq. ( p6[ ) resem- 
bles the grand canonical partition function, and is actually obtained from it 
by the Wick rotation: (^Hi —>■ —iji. 

As an example, let us choose the internal symmetry of the system corre- 
spond to C/(l)gi X • • • X U{l)q^, where the qi are the conserved charges. The 
character of U{l)q- is e"?*'''% so the character of the direct product group is 
exp(i^^^j^ Qi^i). The canonical partition function respecting the exact con- 
servation of charges qi has now the form 



Zq,,...,g^{T,V) = 



1 



.27r i'2-iT 

WJo '''■■■Jo 



X exp 



-« E 

i=l 



^■(r, V",7i,---,7^). 



(47) 



The special cases, Zs and Zb,q,s for a Boltzmannian hadron resonance gas 
are considered in previous sections. 
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6 Summary 

The particle abundances have been computed in the canonical formalism us- 
ing the formulation for the exact conservation of baryon number, strangeness 
and charge in the thermal model of particle production. A good agreement 
with the experimental data of GSI Ni+Ni collisions and of E802 collaboration 
in p — A collisions was reported. 

The good agreement with chemical equilibrium does not mean that the 
particle spectra should follow exactly a Boltzmann distribution since the 
momenta of particles can be severely affected by flow. As an example, a 
model with Bjorken expansion in the longitudinal direction will still have its 
particle ratios determined by Boltzmann factors even though the longitudinal 
distribution is nowhere near a Boltzmann distribution [l8| . 
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